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Abstract. We prove square function estimates in L2 for general operators of the 
form B\D\ + D2B2, where Di are partially elliptic constant coefficient homoge- 
neous first order self-adjoint differential operators with orthogonal ranges, and Bi 
are bounded accretive multiplication operators, extending earlier estimates from 
the Kato square root problem to a wider class of operators. The main novelty is 
that B\ and B2 are not assumed to be related in any way. We show how these 
operators appear naturally from exterior differential systems with boundary data 
in L2. We also prove non-tangential maximal function estimates, where our proof 
needs only off-diagonal decay of resolvents in L2 , unlike earlier proofs which relied 
on interpolation and L p estimates. 



1. Introduction 

In this paper, we generalize the square function estimates from the Kato square 
root problem, to a wider class of operators on L2(R n ; C^), n,N > 1. Previously, 
estimates were known for perturbations of a homogeneous first order constant co- 
efficients self-adjoint partial differential operator D, of the form DB or BD, with 
B being a bounded multiplication operator which is accretive on the range of D. 
Let us first recall how such operators appear in connection with divergence form 
equations. The celebrated Kato square root estimate 



div 
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and B 
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\\\f-&ivAVu\\ 2 « ||Vw|| 2 

for divergence form operators with general bounded accretive coefficients A £ /^(R™; £(C n )), 
were proved in one dimension by Coifman, Mcintosh and Meyer [9] and in full gen- 
erality by Auscher, Hofmann, Lacey, Mcintosh and Tchamitchian [6]. It can be 

written || y/(BD) 2 [u, 0] || 2 « ||RD[u,0]|| 2 , with D = 

acting on vectors of dimension N = 1 + n. This estimate is in turn a consequence of 
square function estimates for the operator BD. This approach to the Kato square 
root estimate was developed in [TJ [8] . 

Operators of the form DB and BD appear not only in connection with divergence 
form operators on R n , but also in connection with divergence form equations on the 
half-space R^_ +n :={(t,x) ; t > 0, x £ R n }, with L 2 (R") or H\R n ) boundary data. 
We recall the following approach to boundary value problems from [3]. Consider 
a divergence form equation 

div t>x A(t, x) V t , x u(t, x)=0, t>0,xe R n , 
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splitting C 



l+n 



C © C n . Write e , ei, . . . , e„ for the 



with coefficients A = a ^, 

c a 

standard basis in R 1+n , with coordinates xo = t, £1, . . . x n , and / x := eo • / for the 
normal component and f n :— f — f ± eo for the tangential part. 

On the one hand, at the level of if 1 (R n ) boundary data w|r«, we consider the 
conormal gradient 

./' 

In terms of /, the divergence form equation is d t f±+div x (c(a~ 1 f ± —a~ 1 bf n )+df n ) = 0. 
The conormal gradient /, with the inward conormal derivative as normal component 
fx, is in one-to-one correspondence with the potential u, modulo curl-freeness and 
constants. Written in terms of /, the curl-free condition is d t f\\ = V x (a~ 1 f ± — a _1 6/n), 
curLj/ii = 0. In vector notation this means that the divergence form equation for u 
is equivalent to the vector valued ordinary differential equation 



7l" 




adtu + bV x u 


h 
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dtf + DBf = 



for / under the constraint / G R(-D) for each t > 0, with D :- 
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div a 




acting 



along R n and B : = 
out to be accretive i 



being a multiplication operator, which turns 



a —art 
ca~ x d — ca~ l b 
and only if A is so. 

On the other hand, at the level of L 2 (R n ) boundary data u|r«, we can write 

fx = (AVu) ± =: div x W| h 

with a tangential vector field v n , for each fixed t, assuming appropriate decay of u at 
infinity, since J Rn f±dx = by the divergence theorem. Inserting this ansatz into the 
divergence form equation and commuting d t and div^ yields div x (d t v n + (AVt,x u )\\) — 
0. Since v n is only defined modulo tangential divergence free vector fields, we may 
choose it so that d t v n + (AVt tX u) n = 0. In vector notation this means that the 
divergence form equation for u is equivalent to the vector valued ordinary differential 
equation 

d t v + BDv = 0, 

for the vector field 

— u 

div" 1 (ad t u + bV x u) 
and d t v + BDv 

Riemann systems. In particular, the n components of v n should be viewed as some 
generalized harmonic conjugate functions. 

Estimates of operators of the form DB or BD are by now well understood, see 
0, [5], Hj. The aim of this paper is to prove fundamental estimates for more general 
operators of the form 

B 1 Di + ,02^2, 

which appear for example when, similar to above, writing a more general exterior 
differential system as a vector valued ordinary differential equation in the variable 
transversal to the boundary. See Section [3j We assume that DiD 2 = but, unlike 
earlier results [1], not that D 2 B 2 BiDi = 0. 



Vx 
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One can view both d t f + DBf 



as generalized Cauchy- 
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We next formulate our results in detail. Consider four operators Di,D 2 ,Bi and 
B 2 acting in the Hilbert space L 2 (R n ; C N ) with norm || • || 2 , where n,N > 1. We 
assume the following. 

• The operators D\ and D 2 are constant coefficient homogeneous first order 
differential operators which are self-adjoint and such that R(-D 2 ) C N(Z?i). 
Assume the partial ellipticity estimates ||A/||2 > ||/||#i(Rn) for all / G 

R(A), * = 1,2. 

• The operators B\ and B 2 are bounded multiplication operators Bi : f(x) i— > 
Bi(x)f{x), x G R n , where B^-) G L oc (R n ; £(C N )), i = 1,2. Assume the 
partial accretivity estimates Re(A/, /) > ||/||| for all / G R(A), i — 1,2. 

Denote by := supj eR (£,^\ r i | arg(£>j/, /)| < 7r/2 the angle of accretivity for £>,; on 
R(A), i = 1, 2. For < a < n/2, define the closed sectors S a+ := {A G C \ {0} ; 
| arg A | < a} U {0} and S a _ := — S Q+ , and bisectors S a := U S a -, as well 
as the corresponding open sectors/bisectors S°±, S 1 ", being the interior of S a ±, S a 
respectively. 

Denote by R(-), N(-) and D(-) the range, null space and domain of an operator. 
Define the operator 

T := B X D X + D 2 B 2 , D(T) := {/ G D(D 1 ) ; B 2 f G D(L> 2 )}. 

The definition of operators ip(tT), <p(tT) in the functional calculus of T, is found 
in Section [21 For a function h on R+ +n , define the (L 2 Whitney averaged) non- 
tangential maximal function 

N.h(x) := sup ( ff \Ks,y)?J^^)'\ * G R", 
<>0 \JJw(t,x) 1 1/17 Is^J I / 

where W(i,a;) denotes a Whitney region around (t,x), for example a?) = 
B(x,t) x (t/2,2t). 

Theorem 1.1. Under the above hypothesis, T is a closed and densely defined opera- 
tor in L 2 (R n ; C^), with spectrum cr(T) C S u , where to := max(w 1 , uj 2 ), and resolvent 
estimates \\{XI -T^W < 1/dist (A, a{T)). 

Moreover, the following estimates of holomorphic functions of the operator T hold. 

• We have square function estimates 

[ rmtT)f(x)\ 2 ^< [ \f(x)fdx, feW), 
Jr» Jo 1 Jr" 

for any holomorphic symbol ip : S° — ± C, u < fi < ir /2, with estimates 
|^(A)| < min(|A| s , |A|~ S ) for some s > 0. 

If furthermore ip\s u+ and i]}\s w - are not identically zero, then the reverse 
square function estimates > hold for all f G R{T) . 

• We have non-tangential maximal function estimates 

[ \N^{tT)f){x)\ 2 dx^ [ \f(x)\ 2 dx, feW), 
Jr» jr» 

for any holomorphic symbol <ft : S° — > C, u < \i < n/2, with estimates 
100) - 1| £ \M S |0(A)| < |A|~ S ; for some s > 0. 
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The estimates in Theorem 11.11 go back to the techniques from the solution of the 
Kato square root problem by Auscher, Hofmann, Lacey, Mcintosh and Tchamitchian [5]. 
The connection between the Kato square root problem and square function esti- 
mates for first order differential operators was developed by Auscher, Mcintosh and 
Nahmod [7]. More directly, both the square function and non-tangential maximal 
function estimates in Theorem 11.11 build on the author's joint work [3] with Auscher 
and Hofmann. 

So far, the Kato techniques have been applied to establish square function esti- 
mates for three main classes of first order differential operators. 

(1) Operators of the form DB and BD, with D being a self-adjoint constant coef- 
ficient homogeneous first order differential operator, and B being a bounded 
multiplication operator which is accretive on R(-D). 

(2) Operators of the form T + BiT*B 2 , with T being a nilpotent (that is T 2 = 0) 
constant coefficient homogeneous first order differential operator, and B±, B 2 
being bounded multiplication operators such that T*B 2 BiT* = = TBiB 2 T, 
which are accretive on R(T*) and R(T) respectively. 

(3) Operators of the form B\D\ + D 2 B 2 , with D\,D 2 being self-adjoint con- 
stant coefficient homogeneous first order differential operators with orthog- 
onal ranges, and By, B 2 being bounded multiplication operators such that 
D 2 B 2 B\D\ = 0, which are accretive on R(-Di) and R(D 2 ) respectively. 

Note that for all these classes of operators, there is essentially only one multipli- 
cation operator B as we are considering a small generalization of the case B 2 = B^ 1 . 
This is in constrast to Theorem 11.11 where the two multiplication operators B\ and 
B 2 are independent of each other. 

The main example of operators of type (T5]) are Hodge-Dirac operators, with F 
being the exterior derivative acting on differential forms, see Section |3j For operators 
of type 02]), square function estimates were proved in [SI Thm. 2.7]. By a simple 
operator theoretic argument, square function estimates for operators of type §T§ 
follow from such estimates of operators of type ([2]), as shown in [SI Thm. 3.1]. In 
fact this argument can be reversed. It was shown in [U Sec. 8.1] that conversely 
square function estimates for operators of type follow from such estimates for 
operators of type (ED). Operators of type ([3]) are nothing but a direct sum of two 
operators of type (CD), and hence square function estimates are immediate as shown 
in [U Sec. 8.2]. The type (JBJ) operators first appeared, in disguise, in the work 
by Auscher, Axelsson and Hofmann [3], where boundary value problems for Dirac 
equations of the form (r + BiT*B 2 )f — 0, Y being the exterior derivative, were 
studied. Similarly to Section [3] here, it was shown in [U Sec. 8.3], that solving for 
the t- derivatives, this equation can be written (d t + {BD\ + D 2 B~ 1 ))Uf = 0, under 
suitable similarity transformation U. It should be noted that in [3] , By and B 2 were 
related in exactly the way so that the associated operator BD\ + D 2 B~ 1 is of type 
([3]) and not of the more general form T considered in Theorem ll.l[ where the two 
multiplication operators are independent. 

Coming to the non-tangential maximal function estimates in Theorem 11.11 these 
build on the estimates by Auscher, Axelsson and Hofmann in j3j Prop. 2.56]. Al- 
though set in the framework with Dirac equations, what was actually proved there 
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was non-tangential maximal function estimates for operators of type (PQ), in the spe- 
cial case when the differential operator is of the form D 



. Even for 
unction es- 



-V 

operators of type (CQ) with more general D, the non-tangential maximal 
timate in Theorem II .11 are new. Also the non-tangential maximal function estimates 
for general operators of type ([2]) and ([3]), which is a special case of Theorem ll.il are 
new. 

The key idea in the proof of the square function estimates in Theorem II. H is to use 
a splitting of L 2 adapted to the operators B 1 D 1 and D 2 B 2 , in which the operator 
B\D\ + D 2 B 2 is triangular due to the assumption D\D 2 = 0. The proof of the 
non-tangential maximal function estimates in Theorem 11.11 is much inspired by the 
proof of [31 Prop. 2.56]. The main difference is that our proof here is a pure L 2 
proof, in that it only requires L 2 off-diagonal decay of resolvents. In [31 Prop. 2.56], 
interpolation theory to prove L p off-diagonal decay, p ps 2, was needed. Another 
novelty is a Caccioppoli type estimate, Lemma HT2"} for operators beyond divergence 
form equations. 

The outline of this paper is as follows. In Section 131 we show how operators of the 
form BiDi + D 2 B 2 arise naturally in connection with exterior differential systems 
systems in R^ +n for differential forms. The special case of one-forms, that is vector 
fields, amounts to divergence form equations. In Section [2] we prove the resolvent 
estimates for the operator T, in Section H] we prove the square function estimates 
for the operator T, and finally in Section [5] we prove the non-tangential maximal 
function estimates for the operator T. The (roadmap to the) proof of Theorem 11.11 
is in Section 

2. Resolvent estimates 

In this section, we establish the basic operator theoretical properties of the oper- 
ator T = BiDi + D 2 B 2 . A fundamental observation for the unperturbed operator 
D\ + D 2 is the orthogonal splitting 

L 2 (R n ;C N ) = (N(A) n N(D 2 )) ffi pi) © = : % © © 
The natural perturbation of this splitting which is adapted to the operator T is 

L 2 (R"; C") = (N(Di) n H(D 2 B 2 )) © R(fii£>i) © R(D^. 
Proposition 2.1. We have a topological (but in general not orthogonal) splitting 

L 2 (R n - C N ) = © B{Hi © n 2 , 

where := (H © H 2 ) D B 2 \U © Hx), B 2 l V := {/ G L 2 ; B 2 f G V} and 
B\V := {BJ ; feV}. 

Proof. We observe that we have two topological splittings 

L 2 = B{K X © {Ho © H 2 ) 

and L 2 = B 2 H 2 © (Ho © "Hi), since B\ is accretive on Hi and B 2 , and hence B 2 , is 
accretive on H 2 . See for example Prop. 3.3]. Taking orthogonal complements in 
the second splitting, we obtain a third topological splitting 

L 2 = B 2 1 (H Q ®H 1 )®H 2 , 

since B 2 l (H © H x ) = (B*H 2 ) L and H 2 = (H © H^. 
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The result is now a consequence of Lemma 12.21 below, with X\ = B{Hi, X 2 = 
U © n 2 , X 3 = B^iUo © U{) and X 4 = U 2 . □ 

Lemma 2.2. Assume that a Banach space X splits topologically in two ways 

X = X 1 © x 2 = x 3 © x 4 , 

into closed subspaces such that X 4 C X 2 . Then X splits topologically into three 
closed subspaces 

x = x 1 © {x 2 n x 3 ) © x A . 

Proof. It is straightforward to verify that these three subspaces are closed and inter- 
sect pair wise only at 0. Also, given x G X, we can write x = xi+x 2 and x 2 = X3+X4 
with Xi G Xi, i = 1, 2, 3, 4. We have x = xi + x 3 + X4, with x 3 = x 2 — £4 G X 3 n X 2 , 
so the three subspaces span X. □ 

Proposition 2.3. For i/ie operator T, the null space is N(T) = T-Lq 2 , the range is 
R(T) = B\R(Di) + R(D 2 ) and the domain is 

D(T) = {u + Ul +u 2 E H% 2 © B{Hi © H 2 ; U\ G Dpi), Ul + u 2 G Dpa^)}- 

Proof. For the null space, we note that / G N(T) if and only if B\D\f = —D 2 B 2 f 
Since B^i n% = {0}, this is equivalent to / G N(Di) n B^N^)- 

Clearly R(T) C -BiR(-Di) + R(D 2 ). For the converse implication, assume that 
fx = B 1 D 1 u 1 G B 1 R(D 1 ) and f 2 = D 2 B 2 u 2 G R{D 2 ) = R{D 2 B 2 ). Write u x = 
u\+u\ E U 2 ® B 2 X (U © Hi) and w 2 = u\ + u| G % © B 2 x {Uq © Then 
{B x D x + D 2 B 2 )(ul + u\) = B 1 D 1 u 1 + D 2 B 2 u 2 = h + f 2) so R(T) = BiR(£>i) + R(D 2 ). 

The result for the domain follows from the facts that T-Lq 2 C D(D\) fl D(D 2 B 2 ) 
and H 2 C D(£>i). □ 

We now express the resolvents of T in terms of the resolvents 

R\ := (I + itBxDx)- 1 and R 2 t := (J + itD 2 B 2 y l 

of B\D\ and D 2 B 2 . It is known that a(B\Di) C S , Wl U {0} with resolvent estimates 
H-R^ll < l/(|£|dist («/£, S^J), and that a(D 2 B 2 ) C S^UjO} with resolvent estimates 
ll-^tll ~ l/(|t|dist (i/t, iS^a)). See for example [SI Prop. 3.3]. 

Proposition 2.4. T/ie operator T is closed and densely defined in L 2 (R n ; C N ). The 
spectrum is contained in the bisector U {0} , u = max(aji, u 2 ), and in the splitting 
L 2 (R n ; C N ) = Hq 2 © B1H1 © T-L 2 , the resolvent has the expression 



(I + itT)- 1 



I ' 

R\ 
(Rl-I)Rl R\ 



i/t $s u \j {o}, 



with estimates + itT) 1 || < l/(|t|dist (i/t, S^)). 



Proof. Consider (/ + itT)u = f with u G D(T), and write u — uq + U\ + u 2 and 
f = fo + fi + f 2 the splitting from Proposition 12.11 Then 

u = fo, ui + itB x D\U\ = /1, 
u 2 + itD 2 B 2 (u 1 + u 2 ) = f 2 . 
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Solving for u, we equivalently have 

«o = /o, v>i = {I + HB 1 D 1 )~ 1 f 1 , 
u 2 = (I + iW 2 B 2 y x f 2 - UD 2 B 2 (I + iW 2 B 2 )-\l + itBxDx)- 1 ^. 

This shows that I+itT is injective with the stated resolvent expression and estimate. 

To show surjectivity, given / G L 2 (K n ; C N ), define u := f + (7 + itD 2 B 2 )- x f 2 + 
(7 + itD 2 B 2 )~ 1 (I + itBiDi) -1 fi. Then reversing the above calculation, shows that 
u G D(T) and (7 + itT)u = f. It follows that 7 + itT is surjective and that T is a 
closed operator. That T is densely defined, follows from the fact that D(T) contains 
the dense subspace D(D 2 B 2 ) n R{D 2 B 2 ). □ 



Proposition 2.5. T/ie adjoint of T is T* 
{feD(D 2 ) ; B* 1 feD(D 1 )}. 

Proof. It suffices to show that if 



= B$D 2 + D X B{, with domain D{T*) := 



(1) D(T) — C : m t— 7- (Tu,v) 

is L 2 continuous, then v G D(_D 2 ) H D(7>i7?*). The splitting for v analogous to 
Proposition 12. II for u, is 

v = v + V! + v 2 G ((H © Hi) n (Biy l (U © %)) © #2% © Hi- 

We need to show v\ G D(D 2 ) and t>i + t> 2 G D(7)i7?*). To this end, let u = u 2 G 
% n D(D 2 B 2 ) C D(T) in (□). Then 

\(D 2 B 2 U 2 ,Vi)\ = 1(^2, U)| < ||«2||2. 

It follows that Ui G D(D 2 ), since (D 2 B 2 )* = B* 2 D 2 . Therefore, for general u G D(T), 
we have 

(Tu, v) = (B 1 D 1 u 1 , v) + (u, B* 2 D 2Vl ), 

so K^DiWi,?;! + v 2 )\ < \\u\\ 2 . Since (B y D y )* = DiB{, it follows that v\ + v 2 e 
D(DiBl). □ 

We end this section with a short discussion of the definition of the functional 
calculus of bisectorial operators. For further details see pQ , where the corresponding 
theory for sectorial operators is readily adapted to bisectorial operators. 

Given a bisectorial operator T in a Hilbert space H, that is a closed and densely 
defined operator T with o~(T) C S u for some uj < 71/ 2 and resolvent bounds 

||(A/-T)- 1 ||<l/dist(A,^), 

there is a natural definition of 4>{T) for any rational function 0(A) which is bounded 
and without poles in S w . Useful such symbols in this paper are for example 1/(1 + 
t 2 X 2 ) and t\/(l + t 2 A 2 ), with scale parameter t > 0. 

If T is not injective, then there is a topological splitting H = Ho © Hi, with 
H = N(T) and Hi = R(T). Indeed, 

7 = (J + itT)" 1 + ztT(7 + ztT)- 1 , 

where the two terms converge strongly to the projections onto Ho and Hi respec- 
tively as t -> 00. Thus T = © Ti and (A7 - T)" 1 = A _1 7 © (A7 - Ti) -1 , where 
Ti := T\ Hl is an injective bisectorial operator in H x . 
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With the Dunford integral 

ij>{T) := mino © I ^(A)(AJ % - T x )-H\ 

the functional calculus is extended to all symbols ip : S° U {0} — > C which are holo- 
morphic on an open bisector S ° D S u \ {0}, with estimates [-0(A) j < min(|A| s , |A|~ S ) 
for some s > 0, so that the integral is convergent in the operator norm on Hi. Here 
the curve 7 = {te ±td ; t G R}, u < 9 < fi, is oriented counter clockwise around S u . 

To obtain 4>(T) as bounded operators on % for general bounded holomorphic 
symbols 0, without decay at and 00, square function estimates 

f°° dt 

J \mtT)f\\kj « 11/11%, /eWi, 

are required. Usually, it suffices to show estimates <, as estimates > for T follow 
from estimates < for T*. Another basic result concerning square function estimates, 
is that if such hold for one symbol ip, then they hold for all ip such that |^(A)| < 
mindAI 5 , \X\~ S ) for some s > 0, and i>\s u ± 7^ 0. 

Given such square function estimates, it follows that il> n {T) are uniformly bounded 
and converges strongly in £(H) whenever sup rigZ+ Ag5 o |^n(A)| < 00 and ip n (\) con- 
verges for each A G S° U {0}. Through such a limiting argument, we construct a 
bounded homomorphism 

taking bounded symbols : S° U {0} — > C which are holomorphic on S° to bounded 
linear operators 4>(T) on %. 

3. Applications to exterior differential systems 

In this section, we show how operators of the form B\D\ + D2B2 appear in con- 
nection with exterior differential systems for differential form. We first fix notation. 
Instead of writing {dxo, dxi, . . . , dx n } for the basis one-forms, we shall keep the nota- 
tion {eo, ei, . . . , e n } from Section [[] for the basis vectors, and we use the terminology 
A;- vector field instead of /c-form, in the euclidean space R 1+n . 

The space of fc-vectors in R 1+n we define to be the ( 1 ^ n ) dimensional complex 
linear space 

A fe R 1+n := span c {e< jl a . . . a e Sk ; < si < s 2 < . . . Sk < n}, k — 2, . . . , n + 1, 

and we let A°R 1+n := C, A 1 R 1+n := C 1+n and A fc R 1+n := {0} if k <£ {0, 1, . . . , n+1}. 

Given a vector v = YTj=o v j e j e A-'-R 1 ^ and a /c-vector w = J2o<s 1 < s k <n w s e s G 
A fc R 1+n , writing s = {si, . . . , Sk} and e s := e Sl a . . . a e Sk , we have in particular the 
exterior product v aw G A fc+1 R 1+n and the (left) interior product ujidG A fc_1 R 1+n 
defined bilinearly using 

j e (j» e {i}us, j£s, JO, j£s, 

e,- a e s := < e,- j e s := < 

[0, j G s, ye{j, s \ {j\)e s \ {j} , j G s, 

where the permutation sign is e(j,s) := (-l)l{»«;J>'iH. Defining inner products on 
A fc R 1+n , k — 0, 1, . . . , n + 1, so that the standard bases above are ON-bases, we have 
that 

v a (•) : A fe R 1+n -> A fc+1 R 1+n and A fc+1 R 1+n -)• A fe R 1+n , 
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are adjoint multiplication operators if v is a vector, that is v G A 1 R 1+ ™, with real 
coefficients. The corresponding differential operators are the exterior derivative op- 
erator 

n 

Vt,* a / := ^ ej a djf, 

j=0 

mapping /c-vector fields / : R 1+n — > A fc R 1+n to k + f-vector fields, and the interior 
derivative operator 

n 

j=o 

mapping k + f-vector fields g : R 1+n — > A fc+1 R 1+n to /c-vector fields. As special cases 
of these operators, we have the gradient and curl, being the exterior derivative acting 
on scalars and vectors (k = and k = 1 respectively), and the divergence being the 
interior derivative acting on vectors. We also note the duality J J (V t , x a f,g)dtdx = 
— //(/, Vt ;X j g)dtdx for compactly supported fields. 

The basic exterior differential system in R 1+n that we want to consider is 

^ J ^t,x -I fk+l = V t ,x A fk-h 

I Vt jX a fk+i = = Vt, x -i /fe-i) 

for a k + f-vector field fk+i and a k — f-vector field fk-i- Two important special 
cases are the following. If k — 0, then the system reads div^/i = = curh j:r /i, 
since /_i = 0. This is nothing but the Laplace equation, written for the gradient as 
in Section [H If k — 1, then the system reads V t)X j f 2 = V t , x fo, V t , x a f 2 — 0. This 
equation is the Stokes' system of linearized hydrostatics, written for the vorticity ji 
and the the pressure fo- 

Consider next a bilipschitz map p : R^ +rt — > f2 C R 1+n , and the system ([2j) in 
fl. We want to pull back this system of equations to R+ +n , and recall therefore 
the following facts from differential geometry. At a fixed point in R+ +n , denote by 
p the Jacobian matrix of all partial derivatives of p. Extend this linear map as a 
A-homomorphism to A fc R 1+Tl , letting 

p(e Sl a ... a e s J := (pe Sl ) a ... a (pe Sk ). 

Given a /c-vector field / : — > A fc R 1+n , we define the pullback of / by p to be the 
/c-vector field 

p*f(t,x) :=p* t ^f(p(t,x))) 

in R+ +n , where p* . is the adjoint of the Jacobian matrix at (t, x). A fundamental 
well known result is that 

(3) V t>x A(p*f)=p*(V t>x Af). 
Less commonly used is the equivalent dual result that 

(4) Vt, x — ' (J P P^g) = JpP' 1 ^^ j g), 
where J p is the Jacobian determinant of p and 

p- l g{t,x) := p {tx) (g{p{t,x))) 
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is the push forward of g by p 1 . Applying ([3]) and (j3J), we find that (T5J) in Q is 
equivalent to 

^ |V M . j {A k+1 (t,x)f k+1 ) = A k (t,x)(V t , x a fk-i), 

[Vt,x a fk+i = = Vt, x J (v4 fc _i(t, x)f k -i) 

in R+ +T \ where /j- := and Aj := J p (p* p*)~ l is the Jacobian determinant times 
the inverse of the metric tensor G = p*p*, extended as a A-homomorphism to A J H 1+n . 

We now show, analogous to the case k = in the introduction, how (jSJ) is 
equivalent to a vector valued ordinary differential equation d t f t + Tf t = 0, for 
general bounded measurable and accretive coefficients Aj(t,x) G £(A J R 1+n ), j = 

— 1, k, k + 1, with an infinitesimal generator T of the form T = B\D\ + D 2 B 2 . We 
use the natural identifications 

A . R n A . + l R n = A fc + l R l + n. . ^ ffi J ^ ^ ^ ^ ^ + = ^ 

Afc -2 R n A fc-l R n = A fc-l R l+n . J fc _ 2 w &Q A + = 

with corresponding splittings of the coefficient matrices so that 

Ak+ifk+i = e a (a k+ if k + h+ifk+i) + ( c fe+i/fe + ^fc+i/fc+i)> 
Afc-i/fc-i = e a (a fe _ij fe _ 2 + 6 fc _i7 fc _i) + ( c fc-i7 fc -2 + 4-i7fe-i), 

and similarly for A^. 

Let ?f A denote the closure of the range of V s a(-) : L 2 (R n ; A i_1 R n ) -»■ ^ 2 (R n ; A J 'R n ), 
and let ^ be the closure of the range of V x j (■) : L 2 (R n ; A J+1 R") -»■ L 2 (R n ; A J R n ). 
Fundamental results are that H 3 A is the null space of V x a (•) : L 2 (R n ; A 3 R n ) — > 
L 2 (R n ; A i+1 R re ), H 3 ] is the nullspace of V x j (•) : L 2 (R n ; A j K n ) -> L 2 (R n ; A- J_1 R n ), 
and we have an orthogonal Hodge splitting 

L 2 (R™;A J R n ) =H{®W J . 

Proposition 3.1. Assume that Aj G L M (R n ; £(A- J - 1 R n © A J 'R n )) are t -independent 
and accretive on W~ l © H 3 A , j = k - 1, k, k + 1. De/ine i/ie A fc " 2 R" © A fe_1 R n © 
A fc R n © A fc+1 R n ^a/wed function f by 

afc-i/fc_ 2 _+ h-ifk-i 

fk-1 

dk+ifk + bk+ifk+i 

fk+l 

Then the exterior differential system (jSJ) /or £/ie A fc_1 R 1+ri © A fc+1 R 1+n valued func- 
tion © / fe+1 zs equivalent to the vector valued ordinary differential equation 



f = 



fk-2 




fk-1 




fk 




fk+1 





dJ+(B 1 D 1 + D 2 B 2 )f = 0, 
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together with the constraint f G R(B 1 D 1 + D 2 B 2 ) for each t, where 
"0 



By 





a k l -a^h 

c k a k x 4 - c k a k l b k 



0' 

Va.j(-) 

-V x a(-) 





D 



2 • = 



V, j (•) 
-V,a(-) 



-1 







-V x aQ 

a k-i -dk-ih-i 

Ck-iO-k-i 4-i — c fc _ 1 a fc _ 1 6fc_i 








V, j (•) 







and 








satisfy the hypothesis in Theorem \l.l[ 
Proof. The equation Wt,x a fk+i = is equivalent to 

dj k+1 - V x a J k = 0, 



V x a/ 



fc+i 



0. 



0. 



The equation Vt, x J (^U-i/fc-i) = is equivalent to 
f V x j K-i7a,- 2 + h-ilk-i) = °> 

\d t (a k _J k _ 2 + & fe -i7 fc -i) + V,j (c fc _i7fc_2 + 4-i7it-i> 
The equation j (^4 fc+ i/fc+i) = A fe (Vt )a; a / fc _i) is equivalent to 

{-V x j (a fc+ i7 fe +_&fc+i7fe+i) = a fc (<9 t 7 fc _i - V^a / fc _ 2 ) + h(V x a 7 fe _i), 
dt(a k +Jk + frfc+iL+i) + V x J_(c fc+ i7fc + 4+i7fe+i) 
= Ck(d t f k -i - Vx A /fc_ 2 ) + 4(V X . A / fc-1 ). 

Written in terms of /, the four evolution equations are 

dtfk-2 + V x . _i (c k _J k _ 2 + 4_i7fe_i) = 0, 

A_i - V x a 7 fc _2 + a* 1 6fc(V s a /jfe_i) + flfe 1 V x — i (afe+i7fe + h+J k +i)_ = °> 
d t fk + V x j (c fc+ i7fc + 4+i7/t+i) - c k (d t f k -i - V x a 7 fc _ 2 ) - 4(V X a 7fc-i) = o, 
K d t fk+i - V x A/*, = 0, 

and the remaining two equations give the constraints V x a f k+ i = = V x . j f k - 2 . We 
next write the evolution equations in terms of /, using f k _ 2 = a k \(f k _ 2 — b k -ifk-i) 
and f k = a k+1 (f k — b k+ if k+ i). The tangential derivatives in the evolution equations 
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D 2 B 2 f. 



that appear with coefficients to the right are 

V x _i (ck-ia k \(fk-2 — bk-ifk-i) + dk-ifk-i) 
-V x a (a^ 1 (/ fc _ 2 - 6jfc_i/fc_i)) 

V x J (CA,-+ia fcH ': i (/yfc — + dk+lfk+l) 

-V x A(a^l 1 (f k -b k+1 f k+1 )) 

The tangential derivatives in the evolution equations that appear with coefficients 
to the left are 





: b k V x a + a fc 1 V x . _i f k 







a k ~\ x jf k ) 



dk^x A /fc_l 

This shows that the evolution equation for / is d t f + (-B1-D1 + D 2 B 2 )f = 0. To show 



SiDi/. 



that the constraint V x a /fc+i 
we note that 



= V x j /fc-2 is equivalent to / 6 R{B x Di + -D2-B2) 



R{B l D 1 + D 2 B 2 ) = B 1 R(D 1 
= B 1 {U k A ®U 



k-l 



R(D 2 ) 
© ((7£* 



7£ 



{H h A +1 © ^ J fe " 2 )) 



L 2 (R n ; A fc_1 R" © ATT 



by Proposition 12.31 and a L 2 Hodge splitting of L 2 (R n ; A fc 1 R n 
Bl 



A fc R n ) adapted to 
□ 



Given Theorem 1 1 . 1 1 and Proposition [3H1 we can proceed as in [SI Thm. 2.3], where 
the case k = was treated, to represent solutions to the exterior differential system 
(IS]) with functional calculus as outlined in Section [21 To this end, define symbols 



-IX 



-IX 



X + (A) := 



-IX 



X"(A) := 



-tx 



Re A > 0, 
Re A < 0, 

Re A > 0, 
Re A < 0, 



t > 0, 



t < 0. 



For t = 0, we obtain bounded spectral projections x ± (T), with x + (T) + x~(T) being 
the projection onto R(T) along N(T). The following result roughly states that the 
spectral subspace x + (T)L 2 := R(x + (^)) is a Hardy type subspace containing traces 
of solutions to (jSJ) in R+ +r \ whereas the spectral subspace x~(T)L 2 : = R(x~(T)) is a 
Hardy type subspace containing traces of solutions to (JSJ) in Rl_ +n , and the operators 
e~ tT x ± {T) are Cauchy integral type operators, giving the value of the function at 
(t, •) from the boundary trace. 

Theorem 3.2. Consider the exterior differential system (jSJ), with bounded, t-independent, 
accretive coefficients A^-i, A k , Ak+i, and the associated operator T = B\D\ + D 2 B 2 
as in Proposition \3 . 1[ Given f$ G x ± (T)L 2 , the function f ~ f defined by 



f±(t,x) := (e~ tT X ± (T))f^(x), (t,x) E R 



l+n 
± ' 



is a solution to (JSJ), with limits lim^o* ||Jt ~ f^lU = an d li m t-»±oo II II 2 = 0. 



Conversely, any solution f ± ~ / ± to (JSJ) R± +n estimates sup t>0 J* 4 



t<±s<2t 



ll/Jl<fa < 
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oo is of this form, and in particular has the stated limits at t = 0, for some 
f G x ± (T)L 2 , and t = oo. 

These solutions have square function, non-tangential maximal function and L^L^ 
estimates 



\\dJt\\ 2 2 tdt * \\N*(f±)\\ 



±i>0 



sup ||/, 

±t>0 



±||2 



II/, 



±112 



The idea of proof is found in |3] Thm. 3.2] and [21 Thm. 8.2]. In particular, 
the estimates follow from Theorem ll.il using the symbol ?p(\) = Ae x (A) for the 
square function estimates, and the symbol 0(A) = e _A x ± (A) for the non-tangential 
maximal function estimates and the L^L^ estimates. We omit the details. 



4. Square function estimates 

In this section, we prove the square function estimates for the operator T in 
Theorem 11.11 We start by simplifying the problem with Lemma 14.11 and we use the 
following operators. 



Pi 

Q] 
P 2 



Q 



(i+t^B.D,: 



\{R t + R l _ t ), 



tB 1 D 1 {I + t 2 (B 1 D 1 ) 2 )- 1 = l(Rl ~ R-t), 



(I + t\D 2 B 2 f)~' = \(R^R\) 

tD 2 B 2 {I + t^D^fr 1 = l 2 {R 2 t ~ R 2 - t )- 

It is known that these operators are uniformly bounded for t > and that square 
function estimates 



WQlfWlj 



\\Q 2 J\\lj< 



fEL 2 (R n ;C N ) 



hold. See for example [3, Thm. 3.4]. 
Lemma 4.1. Let 

Q t := Q t R t Bi. 
Assume that we have square function estimates 



\®tf\\tj< 



for all f G H±. 



Then we have square function estimates J °° ^(t^ 1 )/!!^ < ll/HL / ^ L 2 (R n ; C 
for ip as in Theorem 

Proof. It is known, see [T], that it suffices to prove the square function estimate for 
ip(X) = A/(l + A 2 ). For this ip, we see from Proposition 12.41 that 



^(tT) = §((/+ itTY 1 - (I - itTY 1 ) 












{R 2 _ t -I)Ql + Q 2 Rl Q 



t > 0, 



by writing 



tD 2 B 2 



1 



+ 



tD 2 B 2 



itD 2 B 2 I + itB x D x I - itD 2 B 2 I - itB^ 
itD 2 B 2 tBiDi 
= I — itD 2 B 2 I + t 2 (BiDi) 2 + 



tD 2 B 2 



t 2 (D 2 B 2 ] 



I + itB 1 D 1 
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Since R 2 _ t are uniformly bounded and since we have square function estimates for 
Q\ and Qf, it suffices to prove square function estimates for Q\R\ on BiHi as 
claimed. □ 

We now prove the square function estimates for t using techniques from the 
proof of the Kato square root estimate, following [HI Sec. 4] . 

Definition 4.2. Let T> = IJjgz ^i-i denote the dyadic cubes in R n , with 

V t := {2- j {0, l) n + 2~ j k ; k G Z n }, 2~ j ~ 1 < t < 2~ j ,j G Z. 

Write £(Q) for side length and \Q\ for measure of a cube Q. Given Q G V, write 

MQ) ■= MQ) ■= (2 fe Q) \ (2 k ~ l Q), k > i, 

for the dyadic annuli around Q, where aQ denote the cube with same center as Q 
but with £(aQ) = a£{Q). 

The key tool in the proof of the square function estimates, as well as for the non- 
tangential maximal function estimates, are the following L 2 off-diagonal estimates. 

Proposition 4.3. For any m < oo, there exists C m < oo such that 

\\®tf\\L 2 ( F ) < C m (t/dist (£,F)n|/|| 2 , 

for all f G L 2 (R n ; C^) with supp / C E, and any closed subsets E, F C R n such 
that dist (E, F) := inf{|x - y\ ; x e E,y e F} > 0. 

Proof. These estimates are known to hold for R 3 t , and therefore for R J _ t , P/ and Q J t , 
j = 1, 2, see for example [H Sec. 5]. From this, the estimates for Q t := Ql(I — R\)B\ 
follow as in [HI Lem. 2.26]. □ 

These L 2 off-diagonal estimates enable us to approximate the family of operators 
{@t}t>o by a family of multiplication operators {^t\t>o, where formally j t = 0*1- 
More precisely, we let 

oo 

lt{x)v := ^ ®t(vxA k (Q)){.x), x G Q eV t ,v G C^, 

where XA k (Q) denotes the characteristic function of the dyadic annulus A k {Q). From 
Proposition 14.31 we have the estimate 



oc 



(6) IMk(Q) <^2~ km 2 kn / 2 <C 



k=0 



uniformly for all Q G T> t , t > 0, by choosing m > n/2. 

We also need the following Sobolev-Poincare inequality, see [101 Sec. 7.8]. 



Lemma 4.4. Let 1 < r < n, or r = n = 1, and 1/r* = 1/r — 1/n. Assume that 
1 < q < r* and r < p < oo. Then there exists C < oo such that for all u G H^iTV 1 ) 
and < r < R < oo, we have the estimate 



\u 



u s \\L qi n) < C] Q | V-^ [ ^1 -i 1 1 V ^ | | £p (n) , 



for any convex set Q with diameter R and measure \Q\ and any measurable subset 
S C £1 with measure \S\. 
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Proposition 4.5. We have the estimate 

dt 

|e t /- 7 *£t/||! T < \\nl forallfeHx, 
where E t denotes the dyadic averaging operator 

E t f(x) = E Q f --=^J q f(v)dy, xeQeV t . 

Proof. Let P t ,Qt denote the unpertubed operators P^,Q], that is 
P t := (7 + t 2 D\Y\ Q t := tD 1 (I + t 2 D\y\ 

Write 

®tf ~ ltE t f = 9,(7 -P t )f + {Qt - ltE t )P t f + (%E t )E t (P t - T)f =:/ + //+ ///• 
For the first term we have 

\\ih = \\c8(i-i$)Q t fh<\\Qtfh, 

since (7 + itE>iDi)~ l Bi(t 2 D 2 (I + t 2 D\)) = (7 + itB 1 Di)~ 1 (tB 1 D 1 )Q t , and square 
function estimates for Q t give the desired estimate. 

For the term 777 we note from ([6]) that I^T^ ||l 2 ^l 2 < C. Square function 
estimates for E t (Pt — 7) can be proved as in [HI Prop. 5.7], replacing IT there by the 
operator Di. 

The term 77 we write 

{Qt - lt E t )Ptf = E e *(( P tf - E Q {Ptf))XA k{Q )), on Q E V t . 
Proposition 14.31 and Poincare's inequality in Lemma [4.41 yields 

Qev t \fe=o 

<x> 

£ EE 2 " fcm n p */-^( p */)iiW)) 

oo 

< E 2 - fcm+2fc(n+1) pvpjiii 2(A(Q)) « HtV^/Hl, 

QeD t fc=0 

if we choose m sufficiently large. Since Di is elliptic on R(7J 1 ), the square function 
estimate for 77 follows from that for Q t . □ 

To prove square function estimates for the remaining paraproduct term , ~f t E t f, we 
use the following test functions. For a small fixed parameter e > 0, we define for all 
dyadic cubes Q E D and unit vectors v E C^, the test function 

^:=(7 + (e£(g) J D lJ B 1 ) 2 )- 1 (r/QtO, 

where tjq = 1 on 2Q and suppr^Q C (3Q), ||V77q||oo ^ 1/£(Q)- The parameter e is 
chosen so that the accretivity condition 



Re(v,J ftdx^j > \Q\/2 
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r Q dx 



holds. This is possible since it is known that we have the estimate 

1 

W\jq 

This can be proved by applying [HI Lem. 5.6] with the operator Di, similar to [HI 
Lem. 5.10]. 

Proposition 4.6. We have the estimate 
r°° r drdt 

/ / \E t f{x)\ 2 \ lt {x)\ 2 ^ < 11/111, forallfEL 2 (R n ;C N ). 
Jo Jr» 1 

Proof. By Carleson's embedding theorem, it suffices to show that 

[ m [ \lt(x)\ 2 ^<\Q\, for all Q eV. 
Jo Jq 1 

Following the proof of the Kato square root problem, see for example [HI Sec. 5.3], we 
now do (1) a sufficiently fine sectorial decomposition of C(C N ), run (2) a stopping 
time argument to construct an large sawtooth sub region of the Carleson box Q x 
(0, £(Q)) where the test function fg is paraaccretive, and make (3) a John-Nirenberg 
bootstrapping argument for Carleson measure, to show that it suffices to prove the 
estimate 

f {Q) [ \l t {x)E t r Q \^<\Q\ 
Jo Jq 1 

for all Q G T> and unit vectors v G C^. 

We first note that 

||e t #|| 2 = \\Q 2 t (BJ v Q )-ttQ 2 t R 1 t B 1 (D 1 B 1 )(I+(ei(Q)D 1 B 1 fr 1 (VQv)h 

<\\Qt(Bif v Q )h + t/e£(Q)^/\Q\, 

since Q\ and R] and e£(Q)D 1 B 1 (I+(e£(Q)DiBi) 2 )~ l are uniformly bounded. There- 
fore 



n Q > dt 

J \\®tf v Q \\ 2 j< 11^/5112+ |Q|<|Q|, 



so it suffices to prove 

" m r „,„dxdt 



r dxdt 

/ / \{Qt-itE t )r Q \ 2 ^<\ Q \. 

Jo Jo 1 



Although /q ^ Hi in general, we have that /q — tjqv = —e£(Q)DiBiQl e ^(rjQv) G 
R(-Di). Thus Proposition 14.51 applies to / = fn — tjqV, and it remains to show 

J Q \(& t -l t E t )( VQ v)\ 2 ^<\Q\. 

We have on Q that 

(6t - j t E t )(r]Qv) =J2®t(v{r} Q - l)XA k {Q)) = J2 Qt ( v ^Q ~~ l )XA k {Q)), 

k>0 k>2 

and Proposition 14.31 gives 



IK©* -itE t )( VQ v)\\ L2iQ) < J2m k e(Q)) m (2 nk \Q\) 1/2 < (t/mr\Q\ 1/2 

k=2 
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if m is chosen large enough. We obtain 



r w> r drdt r^> dt 

/ / |(6, - lt E t )( VQ v)\^ < / (t/£(Q)) 2m \Q\T - \Q\ 

Jo JQ 1 JO 1 

and the proof is complete. □ 



5. NON-TANGENTIAL MAXIMAL FUNCTION ESTIMATES 

In this section, we prove the non-tangential maximal function estimates for the 
operator T in Theorem ll.il We first consider the operator D2B2, that is the operator 
T in the special case when D\ = 0. We prove the following result. 

Theorem 5.1. We have non-tangential maximal function estimates 

\\N*WD 2 B 2 )f)\\ 2 <\\f\\ 2 , fen 2 , 

for (j) as in Theorem 

This result was proved in 0, Prop. 2.56] for operators DB with D of the form 

D = ^ q V , which appear in connection with boundary value problems for 

divergence form elliptic system. Below we give a simplified proof for general opera- 
tors of the form DB. Before doing so, we complete the proof of Theorem 11.11 using 
Theorem 15.11 

Proof of Theorem Proposition 12.41 proves that T is closed, densely defined and 
has the stated estimates of spectrum and resolvents. The square function estimates 
for ip(tT) follows from Propositions 14.51 and 14.61 using Lemma [4. 11 Using Proposi- 
tion 12.51 the reverse square function estimates follow by duality. 

To prove the non-tangential maximal function estimates, we first note that esti- 
mates > hold since \\4>(tT)f — /|| 2 — > as t — > and 

supr 1 / U(sT)f\\ 2 2 ds<\\NMtT)f)\\l 
t>o J t 

as proved in [21 Lem. 5.3]. For the estimate <, it suffices to consider the resolvents 
(I + itT)" 1 , since 

~ r 00 dt 

\\NMtT)f-(I + itT)^f)\\l< J ^(^/-(J + ^-Vlllj^H/ll 2 , 

where a proof of the first estimate can be found in [21 Lem. 5.3] and the second 
estimate follows from square function estimates with := 0(A) — (1 -MA)" 1 . By 
Proposition |2.4| it remains to prove 

(7) \W*(Rlf)h< ll/lh, feU 2 , 

(8) \\N*(Rlf)h<\\fh, feB&i, 

(9) ll^(^/)||2<||/lk feBiHi. 
Estimate (I7j) follows from Theorem 15.11 For (JH]), we estimate 

\\N,(R 1 t B 1 v)\\ 2 = \\N,(B 1 (I + iW 1 B 1 )- 1 v)\\ 2 < \\v\\ 2 < \\B lV \\ 2 , v e H Xl 

using that B\ is a bounded multiplication operator and Theorem 15.11 with D 2 B 2 
replaced by D\Bi, for the first estimate, and that Bi is accretive on "Hi for the 
second estimate. 
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To prove @, we introduce the auxiliary non-tangential maximal functions 
N*u(x) ■= sup (V fen r (1+n) If \u(s,y)\ 2 dsdy] , k>0. 

t>0 V J J B(x,2 k t)x(t/2,t) J 

From L 2 off-diagonal estimates for R 2 , as in Proposition 14.31 we get that 

oo 

Nl{R*R]f){x) < ^- km N k M t f){x). 

k=0 

We claim that 

ll^(^/)ll2<2 fcra / 2 ||iVr(^/)|| 2 . 
Thus, choosing m large, we obtain the desired estimates 

ll^ («/)l| 2 <ll^ (^/)l|2<||/|| 2 . 



To prove the claim, assume that N k g(x) > A. It follows that there exists x\ E 
B(x ,2 h t) such that t _ ( 1+ri ) ff B , Xlt i 2 \ x r t i 2t -) \d\ 2 ~ We conclude that for all y G 

B(x\,t/2) we have N®g(y) > A, and therefore 

\{x ; N*g{x) > X}\ < \{x ; M(x {y;5 } Sg{y) > x} )(x) > 2~ kn }\ < 2 kn \{y ; N°g(y) > A}|, 

using the weak L\ boundedness of the Hardy-Littlewood maximal function M , from 
which the claim follows,. □ 

We now turn to the proof of Theorem 15.11 where we write D := D 2 and B := B 2 
to simplify notation. We use the following version of the Caccioppoli estimate for 
equations of the form dtv + BDv = 0. 

Lemma 5.2. Let 77 e C^°(R 1+n ) and consider a function v solving d t v + BDv = 
on a neighbourhood o/supp?]. Then 

i B „( ( ,,)n,( ( ,,)i^<// W( ,,)Piv, (i ,,)ito. 

Proof. Integrating by part twice, using anti-symmetry of d t and symmetry of D, we 
obtain 



(BDv, Dv)r] 2 = - J J (d t v, Dv)r] 2 
= JJ (v, d t Dv)r] 2 + 2 J J (v, Dv)r]d t r) = - jj(v, DBDv)rf + 2 J J (v, Dv)r)d t r) 

= - J J (D(t] 2 v), BDv) + 2 J J {v, Dv)r]d t r]. 

We get the estimate 

Re JJ(BDv,Dv)r] 2 < J J \v\ \Dv\ \q\ \ V??|. 
For each t e R, we have by the accretivity of B that 

\D(rjv)\ 2 dx < Re (BD(r]v), D(rjv))dx. 
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Integrating both sides with respect to t and using the product rule for the derivatives, 
we get 

// \ Dv]W < //(n^xwiv,!)** + // |„| W**- 

Using the absorption inequality, we obtain the stated estimate. □ 

Proof of Theorem \5.1\ To prove the non-tangential maximal function estimates, it 
suffices to estimate the semigroups e~ tDB x + {DB) and e tDB x~{DB), since 

\\N.WDB)f - e' tDB X + {DB)f - e tDB X ~(DB)f)\\ 2 2 

< jf U{tT)f ~ e~ tDB X + (DB)f ~ e WB X -(DB)f\\l^ < \\f\\ 2 , 

where the second estimate follows from square function estimates for DB with 
?/>(A) := 0(A) — e~' A L Moreover, by a limiting argument, we may assume that 

/ G Re- 
consider first e tDB x + {DB)f and a Whitney region W = B(x ,t ) x (t /2,t ). 
Write x + {DB)f = Dv with v G x + {BD)L 2 H D(D). Let P be the orthogonal pro- 
jection onto R(-D), and denote by [Pv] := |-B(x , to) I 1 Jb( Xo t ) P v { x )dx the average 

of Pv G R(-D)- Note that Dv = D(Pv) since P projects along the null space N(D). 
Define the function 

w(t,x) := tijj 1 (tBD)BDv(x) + 4>i(tBD)(Pv — [Pv])(x), (t,x) G R+ +n , 

where 

^x(A) := A-^e"* - (1 - A)/(l + A 2 )) X + (A), 
0x(A) : =(1-A)/(1 + A 2 ). 

In the second term, [Pv] is regarded as a constant function. Since <pi(tBD) is a linear 
combination of resolvents, it extends to a bounded operator Loo(R n ) — > L 2 (B(x , t )) 
as in the estimate (E]). For the first term, the decay of ipi at A = and oo shows 
that we have square function estimates for ipi(tBD). We calculate for t > that 

Dw = De- tBD v = e- tDB x + (DB)f, 
d t w = -BDe~ tBD v, 

using D(j)i(tBD)[Pv] = d t <j)i{tBD)[Pv] = 0. Since d t w + BDw = 0, we get from 
Lemma 15.21 that 

tDB X + (DB)f\ 2 dtdx =l\ \Dw\ 2 dtdx < ff It^wfdtdx, 



|e 

w J Jw J Jw 



with a slightly enlarged Whitney region W = B(x ,2t ) x (t /4, 2t ). Using square 
function estimates, we obtain 

~ r°° dt 

\\N^ x (tBD)BDv)\\l < / \\MtBD)BDv\\l- < \\BDv\\ 2 2 < \\f\\ 2 2 , 

Jo 1 
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so it remains to estimate 4>i(tBD)(Pv — [Pv]). To this end, for fixed to/ 4 < t < 2t , 
use the L 2 off-diagonal estimates for the resolvents of BD to get 

oo 

\\MtBD)(Pv - [Pv])\\l 2[BM)) <J2 2 ~ km W Pv ~ Ml Wo,**0) 

k=0 

oo 

< nk(-m+3n-2n/p+2).n-2n/p+2u v7 ,|2 

~ 11 VV \\L p (B(x ,2k2t ))i 

k=0 

using the Sobolev-Poincare inequality from Lemma H3] with suitable 1 < p < 2 = q, 
Q = B(x , 2 k 2t ) and S = B(x , 2t ). Integrating with respect to t, we get 



t^ n [ I |rVi(*-BD)(Pu - [Pv])\ 2 dtdx 
J Jw 



oo 

<Y,2 k{ ~ m+3n+2) M{\VPv\ p )(x ) 2/p < M{\VPv\ p )(x ) 2/p , 

k=0 



choosing the parameter m large. Therefore, boundedness of the Hardy-Littlewood 
maximal function M on L 2 /p gives 

\\N*(e- WB X + (DB)f)\\l < \\f\\l + ||VP«||1 < 11/111 + \\DPv\\l < 



by the ellipticity of D on R(D). 

The estimate of the semigroup e tDB x (DB) follows from the above estimate upon 
replacing D by —D, since e~ t (~ DB ^x + (~~DB) = e tDB x~(DB). This proves the non- 
tangential maximal function estimate for <p(tDB). □ 
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